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Q ■ 1. Introduction 

o ; 

^SJ , The Pierre Auger Observatory ^ is being used to study cosmic rays with energies larger than 10^^ eV 

with unprecedented precision and statistics. An essential quantity that must be deduced from data 
is the lateral distribution function (LDF) that describes the decreasing of the signals in the water- 
tanks as a function of distance. Knowledge of the LDF is important for the reconstruction of the 
shower core and the shower direction. It can also be compared with model calculations to give useful 
information relating to primary mass. Here we describe how the LDF is measured using the large 
sample of events recorded with the surface detector (SD) array and with a small sample observed 
^ ■ with the fluorescence detectors (FD). For hybrid events, in which SD and FD measurements of the 

, same shower are available, the core position is much better constrained than for SD-only events, thus 

■ providing an important cross-check on the LDF determined from SD measurements alone. 

^ ■ 2. The Fitting Method 

o ■ 



(N 



The water- Cherenkov detectors provide information about Cherenkov photons, which are produced 
when charged particles cross the tanks. The number of Cherenkov photons collected is to a good ap- 
proximation proportional to the energy deposit in the tank. The signal is calibrated in units of vertical 
O I equivalent muons (VEM) jHj. The energy deposit, however, depends strongly on the particle type and 

4-^ ■ the conversion from the Cherenkov signal back to the number of particles in the tank is not obvious. 

^ \ For large tank signals (> 15 VEM) this is not crucial since the uncertainty cr(S') of a signal S (in VEM) 

was determined from data of two detectors positioned 11 meters apart to be (t{S) = 1.06 v^. But 
for small tank signals the number of effective particles, n, is needed since their Poissonian fluctuation 
■ dominates the uncertainty of the signal and is required for the maximum likelihood fit. We have 

H \ introduced a function that gives n for a measured signal S{r): n = P{r, 9, E, A) x S{r) where the con- 

version factor P is called the Poisson factor and is presently assumed to be independent of the primary 
energy, and mass, A, for any distance, r, and zenith angle, 9. The factor reflects the different en- 
ergy deposits of different secondaries and is determined by simulations. Finally we set up a maximum 
likelihood fit to determine the parameters of a trial LDF functional form and, at the same time, the 
position of the shower core, by comparing each tank signal, with its fluctuations, to the value expected 
from the trial function, Sth- L = Hi /p("-i7 Mi) x Tlfc foink^fJ^k) x Hz -P'sat(?^^, Mi) x H™ Fzeroirim, Mm)- 
The individual factors of the likelihood function are determined using the information of tanks at their 
respective distance r. The Poissonian probability density, /p(ni,/Xi), is calculated for small signals 
{Si < 15 VEM). For larger signals the Gaussian approximation is used, /G('^fc,Mfc)■ The effective 
particle number rt of a saturated tank represents a lower limit of the actual signal and we have to 
integrate /g over all possible values larger than n, to estimate the detecting probability of a signal 
larger than n. In case of tanks without a signal we have to sum over all Poissonian probabilities with 
a predicted particle number fii and actual effective particle number < 3. 
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Figure 1. Averaged LDF for sec 6 e [1.2,1.4] (NKG Figure 2. Slope /3 vs zenith angle as function of 
fit and floating slope). energy for SD-only or hybrid geometry. Fits to the 

slope parameter for different energy bins are shown 

3. LDF Measurements 

In contrast to S(IOOO) the shape of the lateral distribution does not change much with energy 
Therefore, the normalisation constant is decoupled from the shape parameter and showers of different 
energies are combined. The LDF was deduced from experimental data using SD-only and Hybrid 
events. Data from January 2004 to April 2005 have been used for the following analysis. High- 
quality events have been selected, which had a successful directional reconstruction with 6 < 60°, 
at least 6 stations with signal above detection threshold, and a core position well inside the SD 
array. Using SD-only events the following LDFs have been investigated: (i) a power law: S{r) = 
5(1000) • (r/1000 m)-'', with a 9 dependent index v = a + b{sec9~ 1), (ii) a NKG-like function 0: 
Sir) = A ■ [(r/rs) ■ (1 + r/r^)]"'' with A = S'(IOOO) • 3.47'', (3 = a + b{sec9- 1) and r^ = 700m (since 
/3 and rg are strongly correlated, we have fixed = 700m and left (3 to vary), and (iii) a function 
used by the Haverah Park experiment 0: S{r) = fcr" if r < 800 m, else {^Y kr~^'i+^'/''^^+^ 
with fixed S, the shape parameter 77 varying with zenith angle, and — 4000 m. These forms were 
fitted to individual events using a maximum likelihood fit of core location and LDF at the same time 
(see section|2l. Two SD-only analyses were performed. First, in a four-parameter fit, besides the core 
location x and y, the slope parameters v, (3 and 77 , respectively, have been varied together with the 
scale factor S'(IOOO). Then a parameterisation of i^, f3 and 77 as function of 9 was determined, which 
was then used in a second analysis fitting only x, y, and S'(IOOO). Figure^shows the averaged LDF 
for secO G [1.2, 1.4] when the NKG assumption is used. For comparison a hybrid derived average LDF 
is shown too (see below for details). An energy dependent threshold effect is apparent at large radii 
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Table 1. Moments (mean, a) of residual distribution of exp. data with various LDFs. Only events with > 6 stations 
were used in the present analysis. For the NKG-like LDF both SD and Hybrid moments are shown, denoted by "sd" 
and "hy" respectively. 
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and reflects upward fluctuations of signals close to the trigger threshold of single tanks. In case of this 
NKG-like function with a free slope parameter, /3, the fit results for E > 10^* eV are given in Table|21 
The energy dependence of f3 is shown in Figure|21and described by (3{E) — a{E) + b{sec9 — 1), with 
a{E) = 2.26 + 0.195 logio(i;/EeV) and b = -0.98. To quantify the quality of the fits the residuals, 
(5 — >S'th)/o'th, and their distributions are computed. For a good LDF the residuals should scatter 
symmetrically around with a = 1. Means and standard deviations of the residual distribution 
are used to compare different LDFs and are given in Table ^ For simplicity only residuals up to 
1500 m are taken into account, resulting in a variance smaller than the expectation value of 1 to avoid 
systematic biases of upward fluctuating signals. The NKG-like function flts the data best, which can 
be seen from the smallest mean residuals and the smallest residual variances. 

Complementarily to the SD analysis, a hybrid LDF analysis was performed. The hybrid reconstruction 
exploits the independent knowledge of the core position to determine the shower axis geometry and 
distance from each detector to the shower core. A maximum likelihood flt (section ISJ is used to 
determine only S(IOOO) and the LDF slope parameter (/?). A NKG-like function was studied in the 
hybrid analysis. 

Hybrid triggers usually include a relatively large number of accidental stations, which in the case of 
low multiplicity events could even outnumber the number of stations that are part of the event making 
the identification of accidentals and candidate stations a difficult task. Therefore, strict quality cuts 
were imposed and only events with at least 6 triggered stations were included in the analysis. That 
reduces the sample size, but at the same time, selecting events with a large number of active stations, 
prevents biasing the LDF slope due to signal fiuctuations. Moreover, as the quality cuts imposed 
on the hybrid analysis are similar to those used on the SD-only 
analysis the comparison is straightforward. Figure [21 shows also 
how (3 varies with 6 for hybrid events. Despite the limited statis- 
tics of the hybrid sample the agreement between SD and Hybrid 
is encouraging. The Hybrid data sub-sample is not large enough 
as to accurately quantify the depencence of /3 with energy at the 
time of writing. 





SD 


Hybrid 


Intercept a 
Slope b 


2.24 ± 0.01 
-0.98 ± 0.02 


2.26 ± 0.17 
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Table 2. Comparison between SD 
and Hybrid analysis on the param- 
eterisation of the LDF slope (NKG- 
like function; E > 10^* eV). 



4. Uncertainty in S(IOOO) 

High statistics are required to accurately describe the LDF, both to reduce the statistical and sys- 
tematic uncertainties. Hybrid measurements (though with much reduced statistics) are a useful tool 
to help identify sources of systematic uncertainty, and it is certain that, as the Auger exposure in- 
creases, the functional form of the LDF will evolve and become increasingly accurate. Increasing 
accuracy will lead to much smaller uncertainties in the reconstructed core position, but the ground 
parameter, ^(lOOO) which is used to determine the energy of the primary CR, is very robust to 
innacuracies in the LDF. It can be shown that by measuring S{r) at 1000 m, fiuctuations in the 
ground parameter due to a lack of knowledge of the LDF are minimised. The result of analysing 
one SD event many times, whilst allowing the slope parameter to vary by ±8% is shown in fig- 
ure 131 This was chosen as a reasonable value for the magnitude of the shower-to-shower fiuctu- 
ations, based on measurements made at Haverah Park where the precision was sufficient to mea- 
sure intrinsic shower-to-shower fluctuations 0|. Analysing the shower with different assumed values 
for the slope parameter, results in a shift in the reconstructed core position, but by choosing to 
measure the ground parameter S{r) at the point where the LDFs intersect, (at 1000 m), the 
effect of the changing slope parameter is minimised. At this point the ground parameter is inde- 
pendent of the LDF. Analysing many showers in this way shows that ropt, the optimum ground 
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parameter has very little dependence on zenith an- 
gle, energy, or the form of the LDF used to recon- 
struct the showers. For example, an analysis of 
~ 500 Auger SD events with energies 10^* '^ eV 
< E < 10^9 gY a^j^jj zenith angles 0° < 6* < 60° 
gives a distribution of Vopt with mean 940 m and a 
rms of 110 m. An analysis of simulated events at 
10^° eV gives a mean of 930 m and a rms deviation 
of 40 m. The prescence of a saturated station (pre- 
dominantly in vertical, high energy events) tends 
to push ropt out by several hundred metres, and af- 
ter an analysis of many showers, at different zenith 
angle and energy, S(IOOO) was chosen as a robust 
ground parameter to measure all showers at. At 
1000 m from the core the mean uncertainty in S{r) 
(across all events) is minimised, and furthermore, 
for the few showers where ropt lies far from 1000 m, 
the uncertainty in S(IOOO) can easily be estimated. 

5. Summary and Outlook 
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Figure 3. LDFs of an Auger SD event, analysed 
with different slope parameters. Red points indicate 
stations with zero signal. The inset shows how the 
optimum ground parameter ropt varies with the array 
spacing. 



The lateral distribution function of EAS observed using the Auger Observatory has been derived. 
Different functions have been tested and it is concluded that an NKG-like LDF describes the data well. 
The dependence of the function on atmospheric depth has been described. It should be emphasized 
that the global shower observables, like the lateral distribution of particles, are not affected by the 
geomagnetic field for zenith angles 9 < 70°. However, for the case of very inclined showers which are 
dominated by muons, the density at ground is rendered quite asymmetric by the geomagnetic field 
and the LDF approach is not longer valid C] . 
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